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Abstract 

An extension of the Standard Model by at least one extra U(l) gauge symmetry has been investigated by many 
authors. In this paper we explore the possibility that this extra U(l) is anomalous. One of the possible signatures of 
this model could be given by the photons produced in the decays of the NLSP into the LSP. 
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1. Introduction 

The start of the LHC has greatly motivated detailed 
phenomenological studies of scenarios which involve physics 
beyond the Standard Model (SM) . Among them D-brane 
constructions in string theory are one of the most promis- 
ing framework in which the SM can be embedded and ex- 
tended. Such brane constructions naturally lead to extra 
anomalous [/(l)'s in the four dimensional low energy the- 
ory and, in turn, to the presence of possible heavy Z' par- 
ticles in the spectrum. These particles should be among 
the early findings of LHC and besides for the above cited 
models they are also a prediction of many other theo- 
retical models of the unification of forces (see [l| for a 
recent review). In [2[ we have considered a minimal ex- 
tension of the Minimal Supersymmetric Standard Model 
(MSSM) with a single extra 11(1)' gauge symmetry in a 
string-inspired setup. We believe that our model encodes 
the key features of the low-energy sector of some of those 
brane construction. In this framework we studied in [3| 
the radiative decay of the next to lightest supersymmetric 
particle (NLSP) into the lightest supersymmetric parti- 
cle (LSP). This kind of process is very interesting since it 
might be the first one where the LSP could be observed 
at LHC [i,|5| and at the upcoming ILC [glHl- 

2. Preliminaries and Lagrangian 

Under suitable assurnptions the LSP in our model 
turns out to be an axino [8| , the fermion component of the 
Stiickelberg supermultiplet related to the anomaly can- 
cellation mechanism (see for details [2|, [sj, |8[). Without 
loss of generality we assume a wino-like NLSP. In the fol- 
lowing we just give the interaction term which involve the 
axino and the wino relevant for our analysis. The interac- 
tion term, written in terms of four components Majorana 



spinorq^, is given by 
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where A2 is the neutral wino, ^5 is the axino, A,^ is the 
photon, 9w the Weinberg angle, go and 32 respectively 
the U{iy and SU{2) coupling constants, A^'^'> the C/(l)'- 
SU{2)~SU{2) anomaly factor and Mz' the Z' mass. The 
rate of the radiative decay (A2 — ^ ^57) is 
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where AM = M2 — Ms, while M2 and Ms are respec- 
tively the wino and axino masses. As we showed in [3|, the 
radiative decay is the most dominant wino decay mode 
with a BR close to 1 (> 94%), so we can use ([2]) to give 
an estimation of the wino mean life time 
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3. LHC Phenomenology 

In order to fall into the WMAP range in the most 
experimentally attractive situation, we considered a light 
LSP (115 GeV < Ms < 150 GeV) and a mass gap of or- 
der am/Ms — 20%, which imply more energetic and 
therefore easier to detect photons. This requirement is 
necessary because the detector resolution increases with 
energy, while at low energy there is an obstruction for 
the detection of photons due to bremsstrahlung, QCD 
background and absorption before the detection from the 
calorimeter iftj. 
Moreover we considered a universal squark mass Mq for 



^The gamma matrices ■y'^ are in the Weyl representation. 
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Figure 1: Number of directly produced wines in function of the 



axino mass Afg and the universal squark mass M 
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the first two squark generations (since under this assump- 
tion they are nearly degenerate) and we assumed flavor 
blindness [lO|. The contribution from the third genera- 
tion squarks is always negligible. In Fig. [T] we summa- 
rize the results obtained in [3| by plotting the number 
of directly produced winos as a function of Ms and Ma 
having assumed 14 TeV of center of mass energy and 100 
fb~^ of integrated luminosity. Since the BR is almost 
close to one this is also the number of photons in the 
final state. The number of photons produced is of the 
order of 10^. In our analysis we follow [9|, ll|- 13[, where 
the NLSP decay in the GMSB framework is controlled 
by the parameter Cgrav ■ If the NLSP lifetime is not too 
long {Cgrav ^ 1) photons originate close to the primary 
interaction vertex ("prompt photons"). In the case of 
large Cgrav and therefore long lived neutralinos the re- 
sulting photons are non-pointing. From now on we fix 
the axino mass Ms — 124 GeV and the universal squark 
mass Mq ~ 3.5 TeV. In our framework the role of Cgrav 
is played by the ratio goA^-^^^Mz'- In the following we 
discuss two different cases: short lived NLSP and long 
lived one. 



3.1. Short life time 

We compare the number of photons produced by ra- 
diative decay with the ones produced by the cascade 
decays of all the other supersymmetric processes. We 
slightly modified the Herwig code 6.5 [IJ] in order to take 
into account the new axino state in the neutral sector. 
It should be stressed that Herwig does not implement 
extra Z' in a supersymmetric framework. This in turn 
implies that the total number of photons can be underes- 
timated due to the lack of sparticles interactions with the 
Z' . However this problem can be overcome by assuming 
a decoupled Z' either because it is very heavy or be- 
cause it is extra-weak coupled. We generated by Herwig 
2-partons— > 2-sparticles events, using about 1 fb~^ of in- 
tegrated luminosity but we have not considered the case 
of SM particles produced directly in the parton-parton 
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Figure 2: Py distribution of photons (in log scale) for 10* susy 
events. 
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Figure 3: Number of generated photons before the preselection cuts. 

interaction. A good discriminant variable of the process 
is the Pt of the photons produced by radiative decay, in 
particular in the region of Pt between 30-80 GeV/c. The 
corresponding distribution is shown in Fig. [2l We denote 
in red the number of 7's radiatively produced from the 
decay of the wino, in blue the number of 7's from all the 
other processes while in black the sum of the two. We 
assumed t^^ — 1.29 • 10~^^s, which is obtainable with 
Mz' ^ 1 TeV and goA^^^ ~ 0.2. We performed the 
same cut on the number of generated photons as in [l3| 
with Pt > 20 GeV and with pseudorapidity I77I < 1.37, 
1.52 < I77I < 2.5, which provides a good way to further 
suppress the SUSY backgrounco. The result obtained 
by using Herwig in generating 10* net events is given in 
Fig. [3l The most important difference between our case 
and the GMSBl sample [9J, |ll|, \l^ is in the number of 
events with zero or two photons in the final state. The 
latter in particular is only 30 in our case. This behav- 
ior can be related to the squark masses we have consid- 
ered. In our case they are about 3.5 TeV, while in the 
GMSBl they are lower than 1 TeV (- 900 GeV). We 
choose the value of 3.5 TeV for the squark masses since 
in this case the number of directly produced winos es- 
sentially depends only on Ms (see Fig. [T]). The number 
of produced squarks is low since they have a high mass. 



^ After having employed the SUSY preselection cut which 
describe later. 



Hence they give a lower contribution to the NLSP pro- 
duction. If we consider hghter squarks, with masses less 
than 1 TeV, there is an increasing in the number of events 
with 2 photons in the final state. In any case the channel 
with one photon in the final state is always the dominant 
one. 

The key point in the analysis is the SM background dis- 
crimination. We considered the same SM background as 
in [13[: events with QCD jets, single gauge bosons (W 
and Z) production, di-Boson, 77 and ti production. In 
order to disentangle the SM background from the signal 
we require a standard preselection cut for SUSY-like sig- 
natures: 

- at least four jets must be present with px > 50 GeV 
{pt > 100 GeV for the leading jet); 

- missing transverse energy i?™'*^ > 60 GeV. 

After having applied these preselection cut we are able to 
reconstruct photons with pT > 20 GeV and \r]\ < 1.37, 
1.52 < I77I < 2.5. The jet-finder algorithm used in our 
analysis is the Durham- type clustering KTCLUS |15| . 
The cut on the ^B™'"^ requires some care. In hadron 
colliders, the initial momentum of the colliding partons 
along the beam axis is not known since the energy of each 
hadron is distributed and constantly exchanged between 
the partons. Hence the total amount of missing energy 
cannot be determined in a straightforward way. However 
the initial energy of particles travelling transverse to the 
beam axis is zero and thus any net momentum in the 
transverse direction denotes missing transverse energy. 
To determine the latter for the i-th not-detected particle 
in each generated event, we considered the following pro- 
cedure: take a transverse direction (perpendicular to the 
z axis) and define the vector 
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whose direction is given by its momentum. Then E™'^'^" 
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^j((£;™'''^)a;)i and analogously for the y component. 
The result is shown in Fig. U) In this plot the num- 
ber of photons from the NLSP decay is that generated 
by Herwig while the number of background photons is 
the reconstructed one [l3|. We note that in the channels 
with N^ > 1 the signal can be easily disentangled by the 
SM background. 

3.2. Long life time 

In the case in which the wino is a long lived particle 
the reconstruction of the emitted photon after its decay- 
ing plays an important role. If the NLSP has a significant 
decay length, a photon will not "point back" to the pri- 
mary interaction point (O in Fig. [SJ but towards a kind 
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Figure 4; Number of generated photons per event for our model 
and for GMSBl and number of reconstructed photons for the SM 
background, after having apphed the cuts described in the main 
text. 
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Figure 5: Schematic diagram of a non-pointing photon, d^^ is the 
distance traveled by A2 before its radiative decay. 2^ is the value 
of the displaced vertex. 



of "virtual" point (O'). The relevant discriminant vari- 
able is Zj, the distance between O and O' . Using simple 
trigonometric relations we have: 



= ^2 - dxy cot ( 



(4) 



with cot 6* = P'^/P^, the ratio between the photon mo- 
mentum along the beam direction (z-axis) and the trans- 
verse momentum. 

By fitting the vertex resolution along z as a function 
of the photon energy for \rj\ < 0.5 (see plot in section 
"Electrons and Photons" in [9|) we derive the functional 
dependence of the standard deviation a for the Gaussian 
distribution associated to z^ on the energy in GeV of all 
the photons as 
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To each z-y we can associate a random number sampled 
from a Gaussian distribution centered on that specific 
value of z-y and with a given by the relation in ([S]) . We 
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considered a wino mean life time r^j ~ 3.463 • 10 
that the distance travelled by this particle before its ra- 
diative decay is always within the ATLAS tracker radius 
(^ 115 cm). In Fig. [5] we show the results of the analysis 



Displaced vertex: \r\\< 0.5 and E > 50 GeV 
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Figure 6: Displaced vertex Zj for photons of B > 50 GeV and 
\ti\ < 0.5: the red hne stands for photons produced by wino decay, 
the black one for photons from the other susy processes. 



for photons of £' > 50 GeV and with \ri\ < 0.5 (photons 
with great angle to the beam direction, because these are 
more easily detected than the ones along the z direction) , 
for an integrated luminosity of 10 fb^^. By considering 
the value of their energy and the contribution of the tail 
of the associated Gaussian distribution, we observe a con- 
sequent enlargement of the distribution of 7's not coming 
from A2, which otherwise would be a spike centered on 
the origin of the x-axis in Fig. [Sj The red line stands 
for photons produced by radiative wino decay, while the 
black one from the others SUSY processes. We can see 
that in the region \zj\ > 3 cm the photons from the ra- 
diative decay are well distinguished from the background 
ones. Since we have assumed AM/Ms — 20%, Ms — 124 
GeV and ta, - 3.463 • lO^i" s, we get from © 



My, 



3.9-10- 
TeV 



(6) 



To satisfy ^ we need either a super-heavy Z' {Mz' — 
385 TeV if goA^"^^ ~ 0.15) or an extra- weakly coupled 
(50^^^' ^ 0.00039 if Mz' ^ 1 TeV) one. This decoupling 
makes our simulation more consistent since in it we have 
neglected the Z' contribution. Cases of extra-weak Z' 
were already studied in the past [l|, |l6| . 

4. Conclusion 

The most important result we get is shown in Fig. [S] 
which gives a distinctive behavior of the photons coming 
from our model with respect to those coming from other 
SUSY processes. The distribution of z^ is centered in 
zero while the half-width is a function of the NLSP de- 
cay length dx^. This means that the radiative produced 
photons are preferably emitted along the NLSP direction 
or more generally at very small angle. Thus the signal 
can be disentangled from the background by applying a 
suitable cut in a region \z^\ < 5 cm. Our result can 
be compared with [9| in which the NLSP is a long lived 
wino-like particle in the particular gauge mediated model 
dubbed GMSB3. The main difference is in the mass gap. 
In our model we considered a mass gap of 20% and an 



LSP of 124 GeV, while in the GMSB3 case the mass gap 
is much higher: the gravitino (the LSP) is almost mass- 
less (m^ = 1.08 • 10"^ GeV) and the NLSP has a mass 
TO^j^ = 118.8 GeV. In this case the energy distribution 
for the photons produced by the NLSP decay is peaked 



on the value of the NLSP mass E 
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while in our case the peak is around 20 GeV. Moreover 
in the GMSB3 ct^, ~ 3.2 m and so the NLSP long lived 
particles tend to escape the detector before their decay. 
This implies a sizable reduction in the number of photons 
detected. 
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